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Abstract 



Using a sample of 3.8 M tp(2S) events accumulated with the BES detector, 
the process tp(2S) — ► ir Jr ir~J/ip is studied. The angular distributions are 
compared with the general decay amplitude analysis of Cahn. We find that 
the dipion system requires some D-wave, as well as S-wave. On the other hand, 
the J/ip - (7r + 7r~) relative angular momentum is consistent with being pure 
S-wave. The decay distributions have been fit to heavy quarkonium models, 
including the Novikov-Shifman model. This model, which is written in terms 
of the parameter k, predicts that D-wave should be present. We determine 
k = 0.183 ± 0.002 ± 0.003 based on the joint m 7Tn - cos6* distribution. The 
fraction of D-wave as a function of found to decrease with increasing 

m-im, in agreement with the model. We have also fit the Mannel-Yan model, 
which is another model that allows D-wave. 

I. INTRODUCTION 

Transitions between bound cc states as well as between bb states provide an excellent 
laboratory for studying heavy quark-antiquark dynamics at short distances. Here we study 
the process ip(2S) — > n+n^J/ip, which is the largest decay mode of the ip(2S) [|l|. The 
dynamics of this process can be investigated using very clean exclusive '0(25') — ► 7t + tt~J/i/j, 
J /ip — ► l + l~ events, where I signifies either e or fi. 

Early investigation of this decay by Mark I |2j found that the 7r + 7r~ mass distribution 
was strongly peaked towards higher mass values, in contrast to what was expected from 
phase space. Further, angular distributions strongly favored S-wave production of TnrJ/if), 
as well as an S-wave decay of the dipion system. 

The challenge of describing the mass spectrum attracted considerable theoretical inter- 
est. Brown and Cahn || and Voloshin M used chiral symmetry arguments and partially 
conserved axial currents (PCAC) to derive a matrix element. Assuming chiral symmetry 
breaking to be small, Brown and Cahn showed the decay amplitude for this process involves 
three parameters, which are the coefficients of three different momentum dependent terms. 
If two of the parameters vanish, then the remaining term would give a peak at high invari- 
ant mass along with the isotropic S-wave behavior. In a more general analysis, Cahn |5j 
calculated the angular distributions in terms of partial wave amplitudes diagonal in orbital 
and spin angular momentum. 

These transitions are thought to occur in a two step process by the emission of two 
gluons followed by hadronization to pion pairs, as indicated in Fig. [1]. Because of the small 
mass difference involved, the gluons are soft and can not be handled by perturbative QCD. 
However, Gottfried || suggested that the gluon emission can be described by a multipole 
expansion with the gluon fields being expanded in a multipole series similar to electromag- 
netic transitions. Including the leading chromo-electric E1E1 transition, T. M. Yan |7j 
determined that one of the terms that Brown and Cahn took to be zero should have a small 
but nonzero value. Voloshin and Zakharov || and later, in a revised analysis, Novikov and 
Shifman || worked out the second step, the pion hadronization matrix element using cur- 
rent algebra, PCAC, and gauge invariance. They were able to derive an amplitude for this 
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process from "first principles." Interestingly, Ref. predicts that while the decay should 
be predominantly S-wave, a small amount of D-wave should be present in the dipion system. 




hadrons 



Hadronization 



V|/(2S) \|f(lS) 

Emission 

FIG. 1. Diagram of ip(2S) — > ir + Tr~J/ip decay process, showing it as a two step process with 
the emission of two gluons followed by hadronization to pion pairs. 

All models predict the spectrum to peak at high mass as it does in ip(2S) — ► 7t + tt~J/iIj 
and T(25) — > 7T7tT(15'). However, T(35) — ► 7r + 7r _ T(lS') [[UJ has a peak at low mass, as 
well as a peak at high mass, that disagrees with these predictions. See Ref. |10| for a list of 
theory papers that attempt to deal with this problem. 

In this paper, we will study the decay distributions of the ip(2S) — > -n + n~ J/ip process 
and use them to test models. The events come from a data sample of 3.8 x 10 6 ip(2S) decays 
taken with the BES detector. 



II. THE BES DETECTOR 

The Beijing Spectrometer, BES, is a conventional cylindrical magnetic detector that is 



coaxial with the BEPC colliding e + e~ beams. It is described in detail in Ref. A four- 
layer central drift chamber (CDC) surrounding the beampipe provides trigger information. 
Outside the CDC, the forty-layer main drift chamber (MDC) provides tracking and energy- 
loss (dE/dx) information on charged tracks over 85% of the total solid angle. The momentum 
resolution is <J p /p = 1.7%\A + P 2 (p m GeV/c), and the dE/dx resolution for hadron tracks 
for this data sample is ~ 9%. An array of 48 scintillation counters surrounding the MDC 
provides measurements of the time-of-flight (TOF) of charged tracks with a resolution of 
~ 450 ps for hadrons. Outside the TOF system, a 12 radiation length lead-gas barrel shower 
counter (BSC), operating in self-quenching streamer mode, measures the energies of electrons 
and photons over 80% of the total solid angle. The energy resolution is oe/E = 22%/ \[E 
[E in GeV). Surrounding the BSC is a solenoidal magnet that provides a 0.4 Tesla magnetic 
field in the central tracking region of the detector. Three double layers of proportional 
chambers instrument the magnet flux return (MUID) and are used to identify muons of 
momentum greater than 0.5 GeV/c. 

III. EVENT SELECTION 

In order to study the process ip(2S) — ► 7i + 7i~J/ip, we use the very clean ip(2S) — > 
7T + 7T~J/if), J/i/j — > sample. The initial event selection is the same as in Ref. fT2|. We 
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require four tracks total with the sum of the charge equal zero. 

A. Pion Selection 

We require a pair of oppositely charged candidate pion tracks with good helix fits that 
satisfy: 

1. | cos^l < 0.75. Here 9 n is the polar angle of the 7r in the laboratory system. 

2. p w < 0.5 GeV/c, where p n is the pion momentum. 

3. pxy n > 0.1 GeV/c, where pxy n is the momentum of the pion transverse to the beam 
direction. This removes tracks that circle in the Main Drift Chamber. 

4. cos 6^ < 0.9, where 9 nn is the laboratory angle between the 7r + and 7r~. This cut is 
used to eliminate contamination from misidentified e + e~ pairs from 7 conversions. 

5. 3.0 < m recoi i < 3.2 GeV/c 2 , where m reco u is the mass recoiling against the dipion 
system. 

6. \ X f /dx \ < 3.0. X f /dx = W^^em^e w here {dE/dx) meas and {dE/dx) exp are 
the measured and expected dE/dx energy losses for pions, respectively, and a is the 
experimental dE/dx resolution. 

B. Lepton Selection 

The lepton tracks must satisfy: 

1. 0.5 < pi < 2.5 GeV/c. Here pi is the three-momentum of the candidate lepton track. 

2. I cos# e | < 0.75, I cos# M | < 0.60. Here 9 e and 9^ are the laboratory polar angles of the 
electron and muon, respectively. This cut ensures that electrons are contained in the 
BSC and muons in the MUID system. 

3. cos# z c ™_ < —0.975. This is the cosine of the angle between the two leptons in the J/ip 
CM, where the leptons are nearly back-to-back. 

4. pi+ or pi- > 1.3 GeV/c or pi+ + p x - > 2.4 GeV/c. This cut selects events consistent 
with J /if) decay, while rejecting background. 

5. For e + e~ candidate pairs: SCE + and SCE_ > 0.6 GeV/c, where SCE is the energy 
deposited in the BSC, or, if one of the tracks goes through a BSC rib or has Pi < 0.8 
GeV/c, the dE/dx information of both tracks in the MDC must be consistent with 
that expected for electrons. The rib region of the BSC is not used because the Monte 
Carlo does not model the energy deposition well in this region. 

6. For pair candidates at least one track must have N hlt > 1, where N htt is the 
number of MUID layers with matched hits and ranges from to 3. If only one track is 
identified in this fashion, then the invariant mass of the pair must also be within 
250 MeV/c 2 of the J/ip mass. 
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C. Additional Criteria 



Fig. 0a shows the m recoi i distribution using the cuts defined above. The shoulder above 
the J/ip peak is caused by low energy pions that undergo it — > \w decay. We impose 
additional selection criteria in order to reduce the amount of mismeasured events from these 
and from other events where the J/ip undergoes final state radiation or where electrons 
radiate much of their energy. These cuts are necessary for comparisons with theoretical 
models. 

1. The 7r's must be consistent with coming from the interaction point. 

2. 3.07 < m recoU < 3.12 GeV/c 2 . 

3. \rrii+i- — mj/ip\ < 0.25 GeV/c 2 , where mi+i- is the invariant mass of the two leptons. 

Fig. ||b shows the m recoi i distribution using all cuts except the additional m recoi i cut (addi- 
tional cut number 2). A total of 22.8 K events remains after all cuts, and the background 
remaining is estimated to be less than 0.3 %. 
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FIG. 2. (a) Number of events versus m reco u, the mass recoiling against the two 7r's, for 
ip(2S) — ► 7T + 7T~ J/if), J/i/j — > events where only the initial selection criteria (see text) are used, 
(b) Number of events versus m reco n where all criteria are used except for the final m reco u cut. 
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IV. MONTE CARLO 

The process is considered to take place via sequential 2-body decays: ip(2S) — > X + J/tp, 
X — > 7r + 7r~, and J/ip — > . The Monte Carlo program assumes: 

1. The mass of the dipion system is empirically given by: 



oc (phase space) x (m M — 4m; 
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2. The orbital angular momentum between the dipion system and the J /if) and between 
the 7r's in the ir + n~ system is 0. 

3. The X and the J /if) are uniformly distributed in cos# in the incoming e + e~ rest frame, 
which is the same as the laboratory frame. 

4. The 7r's are uniformly distributed in cos#*, where 9* is the angle between the J /if) 
direction and the ir + in the X rest frame. 

5. Leptons have a 1+cos 2 9* distribution, where 9* is the angle between the beam direction 
and the positive lepton in the J /if) rest frame. 

6. The J /if) decay has an order a 3 final state radiative correction in the rest frame of the 
J/if). 

A total of 570,000 Monte Carlo events are generated each for the if)(2S) — ► n + 7i~J/if>, 
J/ifj — ► e + e~ and ip{2S) — > -K+n^J/ip, J/tf> — > samples. 

In order to compare with theoretical models, the experimental distributions must be 
corrected for detection efficiency. To determine this correction, Monte Carlo data is run 
through the same analysis program as the data. A bin-by-bin efficiency correction is then 
determined for each distribution of interest using the generated and detected Monte Carlo 



data. This efficiency is then used to correct each bin of the data distributions |L3 |. 

A comparison of some distributions with the Monte Carlo distributions is shown in Fig. |||. 
Fig. [3]a indicates that the m nn distribution agrees qualitatively with the assumed empirical 
distribution ||14j| . Fig. |3|b indicates agreement with the assumed 1 + cos 2 9f distribution for 
leptons in if)(2S) — ► n + 7i^J/if), J/if> — > events. The flat distribution in Fig. ^c is related 
to the assumption that the relative angular momentum between the dipion system and the 
J /if) is zero. However, in Fig. [5]d, which is the cos#*+ distribution, we find a disagreement 
with the Monte Carlo data, indicating that the relative angular momentum of the two 7r's 
is inconsistent with being pure S-wave. 

Fig. [| shows the angle distributions for the l + in the lab; the J /if) in the lab; the ir + 
in the rest frame of the dipion system, and the angle between the normals to the \i\x 
plane and the nn plane. 



arctan 



((X x z) x X) ■ p n -\ 
(X x z)-p^+ 



All distributions are uniform in angle, consistent with the Monte Carlo distributions. 

Since Fig. ||d indicates an inadequacy with the Monte Carlo, it is necessary to correct our 
bin-by-bin efficiency determination in our following studies. We use the Novikov-Shifman 
model (discussed below), which gives a reasonable approximation to the data, to determine 
a weighting for Monte Carlo events so that the proper efficiency is determined as a function 
of cos#* and m W7T . 
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FIG. 3. Various distributions (corrected for detection efficiency) for ip(2S) — » ir + ir~J/'i/j, 
J/vfj — > l + l~ decays, (a) m 1T + 1T - distribution. The distribution is in reasonable agreement with the 
assumed empirical distribution, (b) cos 9f distribution. The assumed distribution is a 1 + cos 2 9f 
distribution. This angle is the angle between the beam direction and the l + in the rest frame of 
the J/tp- ( c ) cos 9x distribution. This is the cosine of the angle of the dipion system with respect 
to the e + e~ direction in the incoming e + e~ CM system. The distribution for Monte Carlo data 
is flat because of the S-wave assumption for the relative angular momentum of the dipion system 
and the J/tp. (d) cos#*+ distribution. This is the cosine of the angle of the tt + with respect to 
the J/tp direction in the dipion rest frame. The Monte Carlo distribution is flat because of the 
assumption that the relative angular momentum of the 7r's is S-wave. The data agree well with 
the Monte Carlo except in (d). 



V. ANGULAR DISTRIBUTIONS AND PARTIAL WAVE ANALYSIS 



In this section, we fit our angular distributions using the general decay amplitude analysis 
of Cahn ||. The ip(2S) and J/ip have J p = 1~ and I GC = , while the dipion system 
has I GC = ++ . At an e + e~ machine, the ip{2S) is produced with polarization transverse to 
the beam. The decay of i/j(2S) can be described by the quantum numbers: 

— * 

I — txix angular momentum 

— * 

L — J/ip X angular momentum 
s — spin of the J/i/j 
s ' — spin of the ip(2S) 
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FIG. 4. Azimuthal angle distributions (corrected for detection efficiency) for 

ip(2S) — > 7r + 7r~ J/ip, J/tp —* l + l~ decays, (a) The eft angle distribution for the l + in the lab. 
(b) The <f) angle distribution for X in the lab. (c) The <p angle distribution for the tt + in the 
dipion rest frame, (d) The distribution of the angle between the normals to the /i/i plane and the 
7T7T plane. 



Defining S = s + I, called the channel spin, then s = S + L = s + I + L. An eigenstate of 
J 2 = s 2 , L 2 , S 2 , and J z may be constructed. Parity conservation and charge conjugation 
invariance require both L and I to be even. 

The decay can be described in terms of partial wave amplitudes, M^l,s, and the partial 
waves can be truncated after a few terms. Considering only M oi, M 2 oi, and M 02 i [P~5|]: 
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The rifTs are measured in their respective rest frames. It is understood that the Mi t L,s are 
functions of m^. The combined 6^ — 9j/^ distribution is given by: 
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A. Fits to ID Angular Distributions 

There are three complex numbers to be obtained. According to Cahn, if the ip(2S) and 
J/ip are regarded as inert, then the usual final state argument gives Mi^s — e * <5 ° ( ' m " r ' r ' ) |^,L,s , |, 
where 5f(m nn ) is the isoscalar phase shift for quantum number I. The phase angles are 



functions of m W7r . If we interpolate the S wave, isoscalar phase shift data found in Ref. JIB 
we find 5q ~ 45°. Also 5® is supposed to be ~ 0. Using these values as input, we obtain the 
combined fit to Eqns. shown in Fig. |5| |TjJ , and the results given in Table Q. Also given 



in Table |I|are the ratios |M 201 |/|M oi| and |M 02 i|/|M oi|. The fit yields a nonzero result for 
I A/201 1 ? indicating that the dipion system contains some D-wave. The amplitude |M 02 i| is 
very small, indicating that the J/ipX angular momentum is consistent with zero. 

Cahn points out that one of the advantages of the process ip(2S) — > 7T + 7T~J/tp is that it 
may allow us to obtain 5q, which is not well measured in this mass range. However we are 
unable to obtain a good fit allowing 5$ as an additional parameter. 

TABLE I. Results of simultaneous x 2 fits to the three 1-D distributions of cos#*, cos#^, 
and cos#* shown in Fig. ||. The phase shifts used are <5§ = 45° and 6% = 0°. The amplitude 
normalizations are arbitrary. Two other fit parameters (not shown) are the normalizations of the 
second and third distributions relative to the first. 



|Mooi| 41.6 ±0.4 ±0.9 

|M 20 i| 7.5 ±1.4 ±1.9 

|M 2i| -0.56 ±0.60 ±0.64 

|M 20 i|/|M 00 i| 0.18 ±0.03 ±0.05 

|Mo2i|/|Mooi| -0.013 ±0.014 ±0.015 

X 2 /DOF 89/111 



9 



- (a) 



-0.6 -0.4 



-0.2 

cos(8*) 



0.2 0.4 



(1» 



-0.8 -0.6 


-0.4 -0.2 0.2 0.4 
cos(9 x ) 


0.6 


O.S 






+ 




(c) 










1 i i t 1 i i i I i i i 1 i t i 1 







-0.8 -0.6 -0.4 -0.2 



0.2 0.4 0.6 0.8 



FIG. 5. Simultaneous fits (x 2 ) to three 1-D histograms of (a) cos#*, (b) cos 9 X , and (c) 
cos#* using Eqns. ||, |l|, and respectively. The phase shifts used are 5q = 45° and 5® = 0°. 



B. Fits to the 2D Distribution 

By integrating Eqn. [| over the <fi angles, we obtain an expression that depends only on 
cos#* and cos 9}^: 
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(5) 



The 2D distribution of cos9^ versus cos9j^ is fit using this equation. We assume 5° = 45° 
and 5 2 = 0°, as was done previously. Using these values, we obtain the fit values shown in 
Table [H]. If we try to obtain <5q, we are unable to get a good fit. 

Fits for different m^, intervals are made assuming 5® = 0° and using values of 5$ that 
depend on the m n7r interval The results are shown in Table (TIT] , along with the values 
of 5q used. The ratios |M 20 i|/|M oi| and |M 02 i|/|M oi| do not show large variations between 
the three intervals, and |M 20 i|/|M oi| is inconsistent with zero for all intervals. 
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TABLE II. Result of the 2-D Likelihood fit to cos 9* versus cos 9* x using Eqn. |5[ The phase 
shifts used are 5® = 45° and <5 2 = 0°. The amplitude normalizations are arbitrary. 



I Mm | 
|M 20 i| 
IM021I 

|M 20 i|/|M oi| 
|M 2i|/|Mooi| 
X 2 /DOF 



13.6 ±0.05 ±0.25 
2.3 ±0.3 ±0.5 
0.05 ±0.16 ±0.19 
0.17 ±0.02 ±0.04 
0.004 ±0.01 ±0.01 
457/437 



TABLE III. Results of the 2-D Likelihood fits to cos 9* versus cos 9* x using Eqn. [| for different 
m^TT intervals. The amplitude normalization is arbitrary. Here the value of Sq used depends on the 
rrin-w interval. 5® = 0. 



m nn Range (GeV/c 2 ) 


0.36 - 0.5 


0.5 - 0.54 


0.54 - 0.6 


Sq used as input 


27° 


42° 


51° 


|Mooi| 


8.36 ±0.04 ±0.23 


7.43 ±0.04 ±0.14 


7.86 ±0.06 ±0.15 


|M 201 | 


1.19 ±0.27 ±0.52 


0.89 ±0.29 ±0.28 


1.37 ±0.37 ±0.56 


|Mo2l| 


0.53 ±0.19 ±0.43 


-0.27 ±0.15 ±0.18 


0.14 ±0.15 ±0.21 


\M 20 x \/\M 001 \ 


0.14 ±0.03 ±0.06 


0.12 ±0.04 ±0.04 


0.17 ±0.05 ±0.07 


M021 \/\M 001 \ 


0.06 ±0.02 ±0.05 


-0.04 ±0.02 ±0.02 


0.02 ±0.02 ±0.03 


X 2 /DOF 


514/437 


608/437 


545/437 


Events 


6186 


7075 


9362 



In comparing the results from Tables ||HTT|, we see that |M 20 i|/|M oi 
0.12 and 0.18 and is at least two sigma from zero. On the other hand, |M 02 i 
between -0.04 and 0.06 and is, in all cases, consistent with zero. 



varies between 
"MqoiI varies 



VI. COMPARISON WITH HEAVY QUARKONIUM MODELS 



A. Novikov-Shifman Model 



A model that predicts some D-wave is the Novikov-Shifman |D| model, which is based 
on the color field multipole expansion to describe the two gluon emission and uses chiral 
symmetry, current algebra, PCAC, and gauge invariance to obtain the matrix element. In 
this model the transition is dominated by E1E1 gluon radiation, so the angular momentum 
of the cc system is not expected to change during the decay and the polarization of the 
ip(2S) should be the same as the J/ip. 
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The model gives the amplitude 



A oc \ q 2 - k(AM) 2 1 + 



2ml 



+ [(AM) 2 - q 2 } 1 



4m 2 
g 2 



cos 2 e - 



(6) 



where q is the four momentum of the dipion system and AM 
parameter k is given by 



M 



M 



The 



(b/6n)a s (fj,)p G (ii) 



(7) 



where b = 9 is the first expansion coefficient of the Gell-Mann-Low function, p G is the gluon 
fraction of the 7r's momentum, which is about 0.4, and k is predicted to be ~ 0.15 to 0.2 
|T9| . From Eqn. |7|, it can be seen that k is expected to be different for ij)(2S) decays and the 
decays of other charmonia, because of the running of a s . The first terms in the amplitude are 
the S-wave contribution, and the last term is the D-wave contribution. Note that parity and 
charge conjugation invariance require that the spin be even. If k is non-zero, it is predicted 
that there should be some D-wave. However, since k is expected to be small, the process 
should be predominantly S-wave. 

The differential cross section is obtained by squaring the amplitude and multiplying by 
phase space. 



dT 



dm W7T d cos 6t 



oc (PS) x A 2 , 



(8) 



where 



PS 



< - 4m 2 )[Mj^ + M; (25) + mU ~ 2(M] ; X, + M 2 (2S) m^ + M^M, 



j/i>' 



4M 2 (2S) 



By integrating over one variable at a time, it is possible to obtain the following ID equations 
for the 171^ invariant mass spectrum and the cos 8* distribution: 



da 

dm^ n 



oc \^{q 2 -4m 2 ) 



+ 0.2k 2 [(AM) 2 - q 2 } 2 1 




, d<T n oc (l.322 - 4.8597k + 5.1577k 2 + 1.18296k fcos 2 01-- 
dcosvt I V 3 



-2.65421 k 2 { cos 2 °* 



+ 0.3738k z ( cos 2 0* 



(9) 



(10) 



The m n7T distribution is fit using Eqn. |^, as shown in Fig. ^. The fit yields k = 0.186 ± 
0.003 with a x 2 /DOF = 55/45. Fitting the cos 9 distribution in the region —0.8 < cos#* < 
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6000 



5000 - 



4000 



3000 - 



2000 - 




1000 



0.4 0.45 0.5 0.55 0.6 

ji + it" Mass (GeV/c 2 ) 

FIG. 6. Fits to the m n7T distribution. The points are the data corrected for efficiency, and 
the curves are the fit results. The smooth curve is the Novikov-Shifman model (Eqn. ^). The 
long-dashed and short-dashed curves are the T. M. Yan model with and without higher order 
corrections, and the dash-dot curve is the Voloshin-Zakarov model (Eqn. |i~3|) . Three of the models 
are nearly indistinguishable. The T. M. Yan model without higher order corrections is slightly 



different. The results are given in Tables IV and VI. 



In! 



500 - 



FIG. 7. Fits to cosO* distribution. The results are given in Tables IV and VI. The points are 



the data corrected for efficiency, and the curve is the fit result using Eqn. 1C . 



0.8 using Eqn. [T^ |2(J, we obtain the results shown in Fig. [7|. The fit yields k = 0.210 ±0.027 
with a x 2 /DOF = 26/40. 

We have also fit the joint cos 9* and distribution (Eqn. ||). This approach does not 
require integrating over one of the variables and is sensitive to any cos 8* - correlation. 
Using this approach, we obtain a k = 0.183 ± 0.002 and a x 2 /DOF = 1618/1482. The 



13 



TABLE IV. 



Distribution 


K 


X 2 /DOF 


m m (Fig. |) 


0.186 ±0.003 ±0.006 


55/45 


cos^ (Fig.0) 


0.210 ±0.027 ±0.056 


26/40 


m n7r vs cos 9* 


0.183 ±0.002 ±0.003 


1618/1482 



results of the different fits are in good agreement and are summarized in Table [TV. 

Using Eqns. || and || where we write Eqn. ^ in terms of S-wave and D-wave parts: 
A = As + Ad, the ratio of the D-wave transition rate to the total rate can be obtained 



R 



D 



J dq 2 J^dcos9(PS)\A D \ 2 
J dq 2 J 1 cos 6 (PS)\ A s + A D \ 



The limits of the q 2 integration are = 4m 2 and g^ aa . = (M^S) — Mj/,^) 2 . For the value 
of k obtained from the joint cos^ — fit, we obtain Rr> = 0.184%. 
The amount of D wave as a function of m n7r has been fit using 

iV(cos^)ocl.0±2^(cos 2 ^-i)±^) (cos 2 £-i) 2 (11) 

The last term corresponds to the amount of D-wave, while the middle term corresponds to 
the interference term ||21|| . The results are shown in Fig. || and in Table [V|. The behavior of 
it as a function of is shown in Fig. |9], along with the prediction of the Novikov-Shifman 
model. 

TABLE V. Fit results to cos#; distribution using a x 2 fit to Eqn. |TI]. The fit also requires a 
normalization term which is not shown. 



Range (GeV/c 2 ) § x 2 /F>OF Events 

0.34 - 0.45 0.319 ± 0.097 ± 0.098 24/37 2016 

0.45 - 0.48 0.085 ± 0.068 ± 0.036 29/40 1995 

0.48 - 0.51 0.144 ±0.045 ±0.033 33/40 3729 

0.51 - 0.54 0.037 ±0.025 ±0.017 35/48 5620 

0.54 - 0.57 0.062 ± 0.022 ± 0.017 44/48 6403 

0.57 - 0.60 0.047 ±0.036 ±0.018 48/48 2959 
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-0.5 0.5 1 -1 -0.5 0.5 1 

cos(O') cos(9*) 
FIG. 8. Fits of cos 9* using Eqn. ^ as a function of m^. The fit results are shown in Table |V|. 
(a) 0.34 < m n7T < 0.45 GeV/c 2 , (b) 0.45 < m n7T < 0.48 GeV/c 2 , (c) 0.48 < < 0.51 GeV/c 2 , 
(d) 0.51 < < 0.54 GeV/c 2 , (e) 0.54 < < 0.57 GeV/c 2 , and (f) 0.57 < < 0.60 
GeV/c 2 . 




0.3 0.35 0.4 0.45 0.5 

m„(MeV/c 2 ) 



FIG. 9. Plot of the interference term, -jt, from Eqn. 11 
prediction of the Novikov-Shifman model for k = 0.183. 



versus m nn . The smooth curve is the 
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B. The T. M. Yan and Voloshin-Zakorov Models 



Other models which describe the m^ invariant mass spectrum are the T. M. Yan Model 
and the Voloshin - Zakarov Model ||. These models are also based on the color-field 
multiple expansion. Yan suggests that the decay can be written as 

(,<T oc (PS) x \(m 2 - 2ml) 2 + — (m 2 -2ml)( m 2 - Ami 

\ J \ 7777 77/ ' Q A \ 7777 77/ \ 7777 77 



dm 




+ 2 ^ 2 U + +0 — , (12) 



where 



K 



M^ 2S) -Mj^ + ml n 

2M^(25) 



The ratio B/A is taken to be a free parameter. The term refers to higher order (HO) 

terms. 

The Voloshin - Zakarov Model calculates the matrix element in the chiral limit, m n = 0, 
and then adds a phenomenological term \m\ 

^-K(PS)xlml-\mlf. (13) 

The m^ invariant mass spectrum has been fit with these models, as shown in Fig. [| As can 
be seen, the Novikov-Shifman and the Voloshin-Zakorov models give nearly identical fits. 
The T. M. Yan model, neglecting higher order terms does not agree as well with the data. 
Including the higher order terms 0, however, gives a fit result which is nearly identical to 
the other two models, as seen in Fig. [| All the results are summarized in Table [V]], along 



with the ip(2S) results from Argus p2| , which used ip(2S) data from Mark II. Argus did not 
fit the T. M. Yan model with the HO corrections, but the the agreement is good for the fits 
they did. 

TABLE VI. Fit Results for the m nn distribution. 



Model 


BES 


Argus - MKII @ 


Novikov- 


K = 0.186 ± 0.003 ± 0.006 


0.194 ±0.010 


Shifman n 


X 2 /DOF = 55/45 


38/24 


T. M. Yan § 


B/A = -0.225 ± 0.004 ± 0.028 
X 2 /DOF = 84/45 


-0.21 ±0.01 


T. M. Yan § 


B/A = -0.336 ± 0.009 ± 0.019 




(HO) 


X 2 /DOF = 60/45 




Voloshin- 


A = 4.35 ±0.06 ±0.17 




Zakorov n 


X 2 /DOF = 69/45 
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C. The T. Mannel - M. L. Yan Model 



Mannel has constructed an effective Lagrangian using chiral symmetry arguments to 
describe the decay of heavy excited S-wave spin-1 quarkonium into a lower S-wave spin-1 



state p4"| . Using total rates, as well as the invariant mass spectrum from Mark II via ARGUS 
|22|] , the parameters of this theory have been obtained. More recently, M. L. Yan et al. |25[ 
have pointed out that this model allows D-wave, like the Novikov-Shifman model. In this 



model, the amplitude can be written ||25 

A oc {q 2 - Cl {q 2 + \q\ 2 ) (l 



2ml 



+ c 2 ml) 



where 



3 
+ 2 



Cl 



Am 2 



cos 2 91-- 



(14) 



— ( 1 + — 



*-> t-l- 1 K, 



-i 



(15) 



and 



1*1 



2m 



V>(25) 



{( m l(2S) ~ ( m ™ + m J/^) 2 )( m l(2S) ~ ( m ™ ~ m J/i>) 2 )} 2 



m„ 



The first term in Eqn. |n| is the S-wave term, and the second is the D-wave term. Note 
that another constant in the effective Lagrangian, has been taken to be zero since it is 
suppressed by the chiral symmetry breaking scale. This amplitude is similar to Eqn. ^| but 
contains an extra term proportional to m 2 . 

We have fit the joint cos#* - m n7T distribution using the amplitude of Eqn. [14| p3| , as 



shown in Fig. 10. We obtain: 



— = -0.49 ±0.06 ±0.13 
9 

— = 0.54 ±0.23 ±0.42 
9 

with a x 2 /DOF = 1632/1481. 

In the chiral limit, g% = 0. If we fit with this value for g 3 , we obtain 

— = -0.347 ± 0.006 ±0.007 
9 

with a x 2 1 DOF = 1632/1482. The results for both cases are given in Table |V11| , along with 
the results from Ref. [|24j which are based on ARGUS-Mark II f2"2fl . The results agree well 
for the (?3 = case. The agreement is not as good for the g% ^ case, but Ref. used 
only the m^ distribution in their fit. In both cases, the x 2 /DOF is large, and there is no 
reason to prefer one fit over the other. 
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i* it" Mass (GeV/c ! ) 
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FIG. 10. Fit of the 2D cos 6* versus m nw distribution to Eqn [L4| (a) The 2D distribution 
projected in m n7T . (b) The 2D distribution projected in cos#*. The points are the data corrected 
for efficiency, and the histogram is the projected fit result. 

TABLE VII. Fit results using Eqn [TJ]. In the second fit, 53 is set to zero. 









91/9 


93/9 


X 2 /DOF 


This Exp. 


-0.49 ±0.06 ±0.13 


0.54 ±0.23 ±0.42 


1632/1481 


Ref. | 


24 


1 


-1.55 ±0.51 


4.07 ± 1.56 


0.87 


This Exp. 


-0.347 ±0.006 ±0.007 





1632/1482 


Ref. [ 


E 


1 


-0.35 ± 0.03 





1.05 




VII. SYSTEMATIC ERRORS 

The systematic errors quoted throughout this paper were determined from the changes 
in the calculated results due to variations in cuts, binning changes in the fitting procedures, 
and changes due to making an additional cut to eliminate background. Cut variations 
include changing the cos 6^ selection from 0.75 to 0.8, changing the pxy n cut from 0.1 to 
0.08 MeV/c, changing the cos^ cut from 0.6 to 0.65, changing the cos# e cut from 0.75 to 
0.7, and changing the m recoi [ cut to 3.05 < m recoi i < 3.14 . 

Fitted results are sensitive to the region of the histogram used in the fitting procedure. 
The changes obtained with reasonable variations in the number of bins used were included 
in the systematic error. 

In addition, the events were fitted kinematically, and a \ 2 cu ^ was made on the fit- 
ted events. Changing the m recoi i cut and cutting on the kinematic fit x 2 determines the 
contribution to the systematic errors due to backgrounds remaining in the event sample. 
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VIII. SUMMARY 



In this paper, we have studied the process ip(2S) — > 7r + 7r~ J/ip. We find reasonable 
agreement with a simple Monte Carlo model except for the distribution of cos#*, which is 
the cosine of the angle of the pion with respect to the J/ip direction in the rest frame of the 
7T7T system. Some D-wave is required in addition to S-wave. 

The angular distributions are compared with the general decay amplitude analysis of 
Cahn. We find that |M 2 oi|/|M oi|, which measures the amount of D-wave of the dipion 
system relative to the amount of S-wave, varies between 0.12 and 0.18 and is at least two 
sigma from zero. On the other hand |M 02 i|/|M 00 i|, which measures the amount of D-wave 
of the J/ip - X system relative to the S-wave, varies between -0.04 and 0.06 and is, in all 
cases, consistent with zero. We are unable to fit for the 7T7T phase shift angle, 5$. 

The distribution has been fit with the Novikov-Shifman, T. M. Yan (with and 

without higher order terms), and Voloshin-Zakorov models. The models give very similar 
fits except for the T. M. Yan model without higher order terms, which gives a poorer fit to 
the data. All fits yield a x 2 1 DOF larger than one. 

In addition, the Novikov-Shifman model, which is written in terms of the parameter k, 
predicts that D-wave should be present if k is non zero. Determinations of k based on the 
cos#* distribution and the joint m nw - cos 6* distribution agree with the value obtained from 
the m n7r distribution. The results agree well with the measurement of Argus using Mark II 
data. However the fit to the joint - cos#* distribution also yields a x 2 1 DOF which is 
larger than one. 

The cos #* distribution has been fit to determine the amount of D-wave divided by the 
amount of S wave, ^, as a function of m^. It is found to decrease with increasing m nn in 
agreement with the prediction of the Novikov-Shifman model. 

Finally, we have fit our m nn - cos 6* distribution using the Mannel-Yan model, which 
also allows D-wave. We find good agreement with their result obtained in the chiral limit 
where g% = using the Mark II data. 

We would like to thank the staff of BEPC accelerator and the IHEP Computing Center 
for their efforts. We also wish to acknowledge useful discussions with R. Cahn, M. Shifman, 
W.S. Hou, S. Pakvasa, Y. Wei, M. L. Yan, and T. L. Zhuang. 
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